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Abstract 

If we want to explain the recently discovered accelerated stage of the universe, one of the option 
we have is to modify the Einstein tensor. The simplest such modification, in agreement with all 
observations, is the positive cosmological constant A. Such a modification will also have its impact 
on local observables and on the propagation of weak gravitational waves. We show here that the 
inclusion of a cosmological constant impedes the detection of a gravitational wave if the latter is 
produced at a distance larger than £crit = (6\/2vr//i/\/5)?^A where ta = 1/\/A and / and h are the 
frequency and the strain of the wave, respectively. £crit is of astrophysical order of magnitude. We 
interpret the result in the sense that the gravitational wave interpretation is only possible if the 
characteristic wave properties are smaller than the non-oscillatory solution due to A. 

PACS numbers: 95.36.+X, 04.30.-w, 04.80.Nii 
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I. INTRODUCTION 



In testing Einstein's theory of gravity, its modifications and ramifications, two important 
sub-areas of research remain to be explored and explained in more detail. The first one 
has to do with cosmology and goes back to the discovery of dark energy ten years ago 
which drives the acceleration of the universe The second one is the possibility to 

detect gravitational waves [3] directly [4] by already operating [b] or forthcoming 6, 7| 
gravitational wave detectors. In order to explain the accelerated universe, we can either 
modify the Einstein's tensor G^i, or try to suitably alter the cosmological energy-momentum 
tensor. The first category encompasses modified gravity theories and theories with the 
inclusion of a positive cosmological constant A js]. This simplest modification is in agreement 
with all observations and, notably, its equation of state p = —p is observationally also 
favored Q|. Once we change the Einstein's tensor to explain cosmological facts, we are also 
forced to accept the fact that the very same parameters will also affect local physics, at 
least in principle. Hence, for instance, the Schwarzschild metric becomes Schwarzschild- 
de Sitter metric where A-effects are also sizable on local scales. Indeed, the theory now 
contains two lengths scales, the small Schwarzschild radius and the large cosmological scale 
ta = ~ Hq^ where Hq^ is the Hubble radius. However, a combination like (rgr^)^/^ 

is of astrophysical order of magnitude and has the meaning to be the distance beyond which 
no bound orbits are possible lOj. Other local effects of the cosmological constant can be 
found in Similarly, the linearized version of the now modified, Einstein's equations will 
include the cosmological constant. These expressions are not any more interpretable as the 
Fierz-Pauli equations [l^ for a spin-2 object. Nonetheless we can still understand them as 
a mathematical approximation for weak fields. Moreover, part of these linearized equations 
will contain the Fierz-Pauli term and therefore the question about gravitational waves in the 
new theory can be also addressed in a meaningful way. What remains to see is how exactly 
the modification of Einstein's tensor influences the propagation of the gravitational waves. 
To this end we solve the linearized equations with A and use them in the energy-momentum 
pseudo-tensor of gravity to study the effect of the cosmological constant. The result, which 
can be formulated in form of a critical distance, is proportional to r\ and depends on the 
frequency and amplitude of the wave and the distance of the source from the detector. 
Although r\ is of cosmological order of magnitude, the small amplitude of the wave arriving 
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on earth renders it possible that the negative contribution of A to the power P is as large as 
the standard oscillatory one. As a result, though the promises to detect gravitational waves 
from most of the systems can be fulfilled, there are some whose gravitational waves detection 
is impeded by A. This is a curious effect of the accelerated universe which gives us also the 
chance to probe the theory of dark energy through gravitational wave detection, provided we 
know the exact distance of the source. The interpretation of our result involves a distinction 
between two tests for gravitational waves: (i) an internal comparison between the two parts 
of the solution (wave and the non-oscillatory part) and (ii) a comparison between the wave 
solution and the cosmological background. 

II. LINEARIZED EINSTEIN'S EQUATION WITH A 



We start with the 



inearized Einstein's equations with A for weak field h^u i-e. the metric 



is Qfiu = Vpty + hfj.!^ where ri^^ is the Minkowski metric (our conventions are like in |l5|): 



= -SttGS^, - Ar/^, (1) 
where we have used the trace-reversed part of the energy-momentum tensor 

SjjLv = Tfj^y — —r]fj_iyT. (2) 
The linearized expression of the Ricci tensor is easily obtained to be 



= ^(DV - d^d^hx, - d^d^hx^ + d^d.h) (3) 



which gives us the linearized equations 

□V - d'^d^hxu - d^d.hx^ + d^d.h = -IQttGS^, - 2Ar]^,. (4) 

This equation is clearly covariant under the local gauge transformation h^^i, h^^ + d^ty + 
dpE^ as imposed by the general diffeomorphic covariance of the Einstein's equations with A. 
Any attempt to make the cosmological constant more dynamical by replacing Krj^i, ^g^iu 
would violate this gauge covariance (in the Appendix we discuss this issue in more detail 
and confirm the validity of (jlj)). This gauge freedom allows us to fix the gauge which we 
choose to be the de Bonder condition: d^h^y = \d^h. The equation to be solved becomes a 
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wave equation with two kinds of inhomogeneities; one the standard source Sfj,u{x), the other 
one a constant term proportional to the cosmological constant: 

□V = - 2Ar]^,. (5) 

Since the equation is hnear we can spht its solution h^^, in two parts 

hf,u = l^Lu + iiiv (6) 
where 7^j, = e^,^(r, u;)e**^"^" + c.c. is the standard retarded solution (written here for a 



monochromatic source at a distance far away from the source [l5j) and ^^j^ solves 0^^,^ = 
—2Krjfj_y. The latter, should satisfy the de Bonder gauge and, in addition, we demand that up 
to a diffeomorphism its asymptotic form is of the de Sitter metric. Both the conditions fix the 
constants a and h and the solution of the homogeneous wave equation QiJ (this is necessary 
to satisfy the de Bonder condition) in the general ansatz where ^^1] = {ar'^+ht^)ri^y. 

In other words, s!ilu is the initial ansatz, supplemented by .^^i"* which guarantees that the 
metric is asymptotically de Sitter, and the de Bonder condition is satisfied. The full solution 



winch i. agree„.e„t with M ( to compare w.th M one ha= to take the g.av.to„ n,a== ™ 
to in [16|) reads, 

2 1 

^00 = -A^^ ioi = gAta;,, ^ij = ^t^Sij + -Ae^^, (7) 

where eij = XiXj for i ^ j and otherwise. These solutions will be used in the energy 
momentum pseudo-tensor t^i, for gravitational waves. 

III. THE ENERGY MOMENTUM TENSOR 



In the absence of the cosmological constant the latter is defined as {0^1, — G^i^)/87rG' |l7 | 
where , again, the index (1) indicates that we expand the tensor in the order 0(h). Taking 
into account that G^u is now modified, the very same procedure can be adopted for theories 
with A leading to 



OTTG 



where tni, is the part defined by 
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Note that t^i, is of the order h^. In agreement with the hnearized equation of motion 
which is of the first order in h and first order in A, the effects of the cosmological constant 
in gravitational waves can be only of the order or Ah. It remains to calculate 

io^ = ^ [-Iho^V^'R^xl + 4? - + Oih^) (10) 

which in the averaged form < ioi > enters the expression for the power of the gravitational 
waves. Making use of the equation of motion Rj}J = — Ar^^^,, we obtain three contributions 
of A to the gravitational Poynting vector ioi, namely 

1 /4 4 2\ .2 1 /10,2 



indicating the different contributions in the same order as in equation ffTOl) . The dots stand 
for oscillatory contributions proportional to 7^ (which average to zero) and the standard 
terms proportional to 7^ surviving the averaging process. The explicit calculation of the 
contribution i?Q- is lengthy albeit straightforward. Assuming the direction of the wave to 
be z, the important quantity for us is < >=< >wavc + < >a where the subscript 
'wave' refers to the standard contribution without the cosmological constant. Taking into 
account that the wave-front moves with the velocity of light (which entitles us to identify 
time with the distance L) one calculates 

^ +03 ^ _^ ^03 ^ _ ^ ^ f03 _ ^ -J r2 /iqN 

SvrG SttG 18 r\ 

where h is either |eii| or |ei2|. Note that due to A, the power 



r 



2-^1 ^ JOi 



dn r 



< r > (13) 



receives a negative contribution. The power is only well defined i.e. positive definite below 
a certain critical distance £crit where the oscillatory character of the solution dominates. To 
calculate this critical distance it suffices to compare the magnitudes of the two contributions 
to < >. The result is 

£„, = 5^4. (14) 

Had we not modified the energy-momentum gravitational pseudo-tensor t^^, to become i^i, 
in equation ([8]) the contribution of A would be bigger and the critical distance smaller by 
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a factor 0.8 which would increase its relevance. In interpreting the above result we empha- 
size that there is little doubt that a modification of the Einstein's tensor will change the 
linearized version of the Einstein's equation (in our case with A this is equation (jlj)). One 
could also paraphrase th.s in saying that the Newtorrian Limit will ehange R. As a con- 
sequence, the solution will now contain an oscillatory and a new contribution originating in 
the modifications of Gfj^i, (proportional to A in our case). The interpretation of gravitational 
waves as ripples on spacetimes can be only maintained if the oscillatory solution is more 
sizable than the non-oscillatory one proportional to A. The result in (fT^ refiects exactly 
this fact. One can also view this result from a more formal perspective. Even though the 
cosmological constant is not part of the energy-momentum tensor, one can nevertheless, 
formally, absorb it there as evident in (jlj). It is then obvious that A will be always a source 
for the metric, gravitational waves to which it contributes, not excluded. 

Notice that what we are really comparing is the averaged solution proportional A with 
the averaged wave component of the solution. We then say that the wave character of the 
solution is lost when both are comparable. 



TABLE I: Sources of gravitational waves for LIGO from references [2^]. AIC means accretion 
induced collapse. For P we have used geometrized units G = c = 1. 



System 


f [Hz] 


h 


Distance 
[pc] 


■^crit 

[pc] 


dP/dn 
A = 


dP/dn 
A/0 


NS/NS binary 


100 


1 X 10-23 


10^ 


12.9 X 10^ 






BH/BH binary 


100 


1 X 10-22 


2 X 10^ 


12.8 X 10^ 






Collapse and explosion 
of Supernova 


20 


4.1 X 10-23 


10^ 


10.6 X 10^ 


1.11 X 10-12 


1.18 X 10-13 


NS formed 
from AIC 


450 


8 X 10-23 


10^ 


46 X lO'^ 


2.15 X lO-'^ 


2.05 X 10-^ 


NS/NS binary 


1000 


1 X 10-20 


2.3 X lO'^ 


12.8 X 10^° 


8.79 X 10-^ 


8.79 X 10-"^ 


Stellar collapse 
Centrifugal hang up 


100 


1 X 10-22 


15 X 10^ 


12.9 X lO'^ 


3.74 X 10-1° 


3.69 X 10-1° 
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IV. PHENOMENOLOGICAL RESULTS 



In exploiting (fT^ phenomenologically, we first point out that the gravitational waves 
arriving on earth are indeed weak which is exactly the reason making their detection difficult. 
They can also be considered weak over the largest part of the distance they travel to earth. 
Therefore, even if (IT^ is an approximation, it is a rather good one. Secondly, the very same 
fact that the waves arriving are weak makes Cent of astrophysical order of magnitude in spite 
of the large value of r\. To see that, let us take some typical values: / = 0.38 x 10~^Hz 
and /i = 40 X 10~^^. We obtain £crit = 1957pc. The values taken are for the white dwarf 
binary system WD 0957-666 whose distance from earth is expected to be roughly lOOpc. 
In this case even though the critical distance is of astrophysical order of magnitude, the 
gravitational waves from the white dwarf devil's system will be seen as its distance from 
earth is smaller that the critical one. The detector sensitive to the values of frequency 
and amplitude (strain) would be in this case the forthcomiiig space-located LISA detector 
j6|. Another example is the collapse of rotating stare cores 19 1 suited e.g. for the planned 



Euro detector 



20[ |. With the characteristic amplitude hf. ~ 10 the frequency / 



900IIz and the relation = ^J^xfrhj2 2l| with r 10~'^s the duration time, we obtain 
-^crit = 1.16Mpc. The maximally allowed distance from earth is supposed to be c? = 15Mpc, 
which implies that the range for the gravitational waves to be detected is much smaller 
than d if A enters the Einstein's equations with the value needed to explain the accelerated 
universe. In tables I and II we have listed three kinds of examples for the LIGO and LISA 
detectors, respectively. In examples where the wave is not monocromatic, we pick up one 
frequency and the corresponding amplitude. The first two entries serve the purpose to 
demonstrate that indeed according to f|T^ the detection of some gravitational waves will be 
impeded by A. The next two examples show that the two contributions to the gravitational 
Poynting vector can be of the same order of magnitude reducing thereby the power of the 
gravitational wave. Finally, and this is the majority of cases, the last two examples show 
that the effect of A can be also negligible. This allows us to conclude that constructing 
a more exhaustive map of all sources whose gravitational waves will not be seen on earth, 
provided the cosmological constant is the right explanation of the accelerated universe, is 
a worthwhile undertaking. Maybe in the near future we will enjoy to see the connection 
between dark energy and gravitational waves which is not only important for the latter, but 
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converts the gravitational wave detectors partly also in experimental devices to check dark 
energy. A better knowledge of the distance of the source is here a crucial ingredient and 
would require an improvement. Indeed, the lOOpc which appears so often in table II seems 
to be only an order of magnitude estimate. If its estimate goes up by a factor 2 — 5, several 
sources might fall into the category whose gravitational waves will not be see due to A. Thus 
the good knowledge of characteristics of the source are of utmost importance for the critical 
distance. 



V. INTERPRETATION 

In this section we argue that due to the appearance of A in the Einstein tensor (and not 
in the energy-momentum tensor) two tests of gravitational waves are required. The first one 
(the cosmological test) is more standard and is due to the interpretation of the gravitational 
waves as ripples on spacetime. Here A appears in the solution of the cosmological back- 
ground. The other test based on (HM is between two solutions: the oscillatory part versus 
the non-oscillatory proportional A. It appears at the first glance that this is the same, es- 
pecially as at present epoch our universe is dominated by A. To see that these two tests are 
different, imagine a universe with a non-zero cosmological constant (say, of the same value 
as in our universe) where, however, the cosmology is dominated by the background density 
(i.e. we can neglect A is the cosmological equations). The crucial point is now that the 
second test relying on (IMI) would be still required and its outcome would be just the same 
as presented in the tables. The effect of Dark Energy models which modify the cosmological 
fluid equation (i.e. the ingredients of the model enter the energy-momentum tensor and not 
the Einstein tensor) could be probed only in a cosmological test which is not what we have 
done here. This is also true for effects which rely on general equation of state distinguishing 



the Dark Energy models 221]. Such a distinction does not make a difference if the model 



modifies the Einstein tensor or the energy-momentum one. 
A. A in cosmology and local physics 

We seem to be biased by the the name "cosmological constant" which has instilled in 
some the impression that A is good for cosmology and nothing else. As mentioned in the 
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Introduction, this is not correct, but it makes sense to look at it from a different perspective. 
Consider, for instance, a gravity theory defined by the more general action 

5 = K y d^x^f{R) (15) 

where f{R) is more complicated than the standard f{R) = R- Any new constant which 
appears in f{R) will also appear in the calculation of any effect of local physics. It would 
be hard to argue that no other constant apart from the Newtonian one (Gn) can enter 
local physics if the gravitational theory is described by f|T5|) . Equally, it would be hard to 
argue that no constant other than G^ can affect the solution of the linearized versions or 
that the effects due to the new constant are coordinate effects. In view of f lT^ we would 
be forced to interpret the gravitational waves anew since part of the solution would involve 
the new constants entering f{R). The situation with A is just a special case of what we 
outlined above. What is really required in cases where the standard G^,^ is modified, is to 
pay attention to the interpretation of gravitational waves. 

First let us note that no local effect of A is per se a coordinate effect. One cannot get 
rid of the cosmological constant A in deriving effects on local physics, in general, and in the 
linearized version of Einstein's equation, in particular. A is an integral part of the modified 
Einstein's tensor G^u and not of the energy-momentum tensor. This implies that A will 
appear locally, in principle, everywhere where gravitational effects are considered. More 
specifically, we have two equations. The first, 

G,. = R,. - \g,uR + f^g,u = kT^~ (16) 

for the universe defines the cosmological background and the Hubble flow and the second 
one where A affects the expansion of the universe, and 

= R,, - \g,uR + lyg,u = ^tJ^^' (17) 

for the local physics including the gravitational waves. In spite of the fact that A appears 
already in the cosmological part f|T6|) . it makes its entry once again in calculating local 
effects. In a concrete example, a star is part of the Hubble flow (expansion) where A already 
contributes, nevertheless the Schwarzschild-de Sitter metric (see e.g. 3]) will again contain 
terms with A which is a local effect of this constant. Another way to see it, is the Newtonian 

nn 

Limit. A survives the Newtonian Limit (|18l. l25l|) 



r, 1 r 



2 



$ = - Vs = GnM, rA = 1/VA (18) 
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because locally its effects are not coordinate effects. Note that the background density ph in 
^.Universe -^^ gquatiou ^ doCS UOt appear in the Newtonian Limit nor in the Schwarzschild- 
de Sitter metric. The Schwarzschild-de Sitter metric is constructed with the boundary 
conditions that its asymptotic form de Sitter. For the latter we could also demand the 
asymptotic form to be the metric of the cosmological background, i.e. Friedmann-Robertson- 
Walker metric 26]. In such a case, A would eventually enter twice, once in the differential 
equations through f|T7|) and the second time through the boundary condition where A is part 
of cosmological background metric. This clearly shows its double role, due to the fact that 
it is part of Gfj,^. 

Hence strictly speaking we have to compare the oscillatory part of the solution with the 
non-oscillatory one ([7]). This could be done in a simplistic way taking the amplitude of the 
oscillatory part h and comparing it with L'^/r\ which comes from ([7]). This comparison 
would conceptually not be very different from what we obtained in f|T4|) . However, this 
way the critical distance would come out smaller than ( fT4l) (of the order ^/hr/<^) as f|T^ is 
suppressed in addition by fr/<^. It makes therefore more sense to make a more sophisticated 
analysis as done in section 2. If the distance of origin of the gravitational wave is larger 
than £crit; the non-oscillatory background is larger than the actual wave and therefore the 
wave interpretation untenable. 



In 271] a condition for the validity of the picture of a gravitational wave as a ripple 
on spacetime is given. Essentially it states that the wave length must be much smaller 
then the curvature background TZ. Since the part of solution which is connected to A is 
non-oscillatory we cannot make such a direct test. In case of a non-zero graviton mass m 
trigonometric functions cos{mt), sin(mt) would enter the solution ([7]) as shown in [161]. Then, 
the above criteria would apply. In our case, we could compare the change of by calculating 
6C,fj,uL/TZ ~ 0{L^/r\) which is much smaller than one as long as L is of astrophysical order 
of magnitude. The reason why the oscillatory solution becomes comparable to the non- 
oscillatory is because the amplitude of the oscillatory is small. 

B. LOCAL VERSUS GLOBAL TEST 

The central point of our interpretation is that given any modified Einstein tensor, there 
has to be two tests of the gravitational waves: 
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1. Global cosmological test: This test can be presented without any recourse to the details 
of the modified Einstein's tensor. The test consists in global comparison between a 
given background cosmology and the wave solution /i™^^'' of the gravitational wave. It 
is crucial to realize that we are comparing here only part of the full solution of the 
gravitational wave, namely the wave part (in the case of A we called it 7^;/). 

This pre-assumes, however, that the wave part is the dominant part of the full local 
solution h^i, = + h^^^^ where the Ii^^^^ is due to the modifications in the Einstein's 
tensor (^^j, in our specific case). Such a comparison of /i™^™ with the cosmological 
background is exactly the Misner-Thorne- Wheeler method |27| mentioned above. One 
does not compare h^^, or h^^^^ with the cosmological background, but only the wave 
part hj[^'^'^ such that all quantities needed for such a test like wavelength, background 
curvature etc are well defined. The cosmological constant enters here only through 
cosmology i.e equation [161 In the present paper we are not concerned about this global 
cosmological test. We rather assume that all sources for LIGO and LISA passed this 
test already. 

There will be also other effects whose root can be traced back to the cosmological 
background. One of them will be the direct appearance of such a background in the 
propagation of gravitational waves for large distances. This effect can be taken into 



account by expanding the Einstein's equations around the de Sitter metric [28|| . Such 
a procedure to include cosmological effects is not general (indeed a general procedure 
does not exist), but is for the present epoch of the universe which is dominated by 
A. The most important effect is the exponentially decay of the wave j29|. As long as 
r ^ ta (or T ^ Ta = ta) we can, however, still rely on the expansion around the 
Minkowski metric. 

2. Local test: In the global test above we made the assumption that /ij^^™ is the dominant 
part of the full solution. This has to be done in a more quantitative manner i.e. we 
have to device a local test which will decide when the wave character is dominant. In 
this local test A enter through equation (17). Our suggestion for such a test is based 
on the power P as explained in the text above. This test is rather conservative as 
other, more 'naive' tests yield a smaller £crit- 

For a better understanding the difference of the two tests, let us visualize a universe 
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(or, equivalently, an earlier epoch of our universe) where the cosmology is dominated 
by the background density and not A. We could then drop A in (16), but not in (17). 
The appearance and relevance of A in the local test would pertain i.e. the results of 
such a test would be the same in any epoch of the universe or any universe with the 
same A (and different background density). 



TABLE II: Sources of gravitational waves for LISA from references 2j]. The last entry is a special 
white dwarf binary. For P we have used geometrized units G = c = 1. 



System 


f [Hz] 


h 


Distance 
[pc] 


-^crit 
r ^ 

[pc] 


dP/dn 
A = 


dP/dn 
A / 


X-ray pulsar binary 
4U1626-67 


7.9 X 10-"^ 


6 X 10-24 


8000 


61 






X-ray pulsar binary 
4U1820-30 


3 X 10-3 


2 X 10-23 


8100 


773 






White dwarf binary 
WD 2331-F290 


1.4 X 10"^ 


2 X 10^22 


100 


360 


1.3 X 10-3^ 


1.2 X 10-31 


White dwarf binary 
WD 1101-F364 


1.6 X 10""^ 


2 X 10-22 


100 


412 


1.7 X 10-31 


1.6 X 10-31 


White dwarf-B star 
KPD 1930+2752 


2.4 X lO"'' 


1 X 10-21 


100 


3090 


9.57 X 10-30 


9.56 X 10-30 


RXJOSO 
6.3+1527 


6.2 X 10^3 


4 X 10-22 


300 


32 X 10^ 


9.2 X 10-2^^ 


9.2 X 10-2^^ 



VI. CONCLUSIONS 

Any gravity theories with modified Einstein's tensor will require a re-interpretation of the 
picture of gravitational waves as ripples on spacetime. The solution of the linearized new 
Einstein's equations will contain the oscillatory part (wave) plus a new component associated 
with the extension of Einstein's tensor. The latter will not be oscillatory, in general. It is 
clear that a suitable comparison between these two solutions is due in order to be able to 
say when the wave picture can be maintained. In this paper we suggested such a comparison 
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by using the gravitational energy momentum tensor associated with observables. Applying 
the method for a theory with the cosmological constant, we deduced a maximal distance 
beyond which the wave picture loses its meaning. This makes a direct connection between 
gravitational waves and theories with A explaining the accelerated universe. 



Appendix: The Veltman Lagrangian 

It is instructive to re-derive the same linearized equations as in (jlj) and to cast a brief 
glance at the .ea.on why the te.. p„po.tio„a, A . not dynanncg, i.e. ptopo.t.onal to 

For this purpose we evoke the Lagrangian given by Veltman in |30|] which reads 

Ch = -2A + - Ih^/sh'"' + hhj - ^d.Kpd^h'''' + ^d^hd^h - ^^dphd^hf^' (A-1) 

The part proportional A is not gauge invariant under h^j_y h^y + d^e^ + d^efj,. Indeed, one 
obtains under this transformation 

2A (^1 + - ^ V^"^ + ^hh 

2A (l + - Ih^/sh'^'' + hh + dPe, - K^d^^e^ + \hd''e}j (A-2) 
The formal condition for the gauge invariance to hold up to total derivative is obviously 

KpO'^e" = hidPe, (A-3) 
The correct gauge invariant Lagrangian is simply 

= Ch + 2A{l/8hh - l/Ahaph""^) (A-4) 
In vacuum, the Euler-Lagrange equations according Ch come out to be 

Dh^"" - ri^^'Dh + ri'^^dM'"" + d^d^h - d^d^h'"' - d^d'^h"'' = 2Ar/^'^ - 2A (^h>"' - ^ri>"'h 

(A-5) 



The second term on the right hand side of ( lA-Sp is due to the the non-gauge invariant terms 



in the Lagrangian. Dropping this term results in equation of motion which we had before 
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i.e. @. (this is equivalent to use the gauge invariant Lagrangian (1A-4P . Indeed, taking the 
trace of 



Dh'"' - r]^"'Dh + r]f"'dM'"" + d^'d'^h - d^d^'h'"' - d^d'^h''^ = 2Ar]^"' (A-6) 

and multiplying with 77^,^ we can replace —r]^i,Dh = —rj^ydadah"^ + 'iKr]^^^ in flA-6P to arrive 
at (jl]) in vacuum ( with matter the steps to obtain (jlj) are similar). This shows once again 
that equation (jll) is correct. 

In passing let us make a comment regarding the gauge invariance of (lA-ll) . Taking the 
divergence of equation (1A-5P gives us 



d^K. = ^d,h (A-7) 

which is actually the de Bonder condition, now not as a gauge fixing, but as a result of the 
equation of motion (this is in analogy to the free massive vector case where in spite of 
the loss of gauge invariance the equation of motion gives us the Lorentz gauge df^A^^ = 0). If 
we put this equation into the Lagrangian ( lA-ll) . then according to ( ]A-2I) and ( ]A-3[) the total 



Lagrangian would come out now gauge invariant up to total derivatives. Obviously, this is in 
contradiction with our previous result and the resolution of this seemingly different results 
is that it is not permissible to use equations of motion (or a part of them) in the Lagrangian 
itself. Similarly, we cannot use a gauge fixing in (1A-4P without changing physical results. 
For instance, if we use the traceless gauge h = in flA-4p . the term {1/A)dyhd'^h in flA-4p 



would be absent. Such a term under variation of the action gives {l/2)rifj_^nh which turns 
out to be crucial to obtain the equation (jll) as explained above. 
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